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Abstract 



We study the nature of the 3 He superfluid transition from the normal to the 
planar phase, which is expected to be stabilized by the dipolar interactions. 
We determine the RG flow of the corresponding Landau-Ginzburg- Wilson 
theory by exploiting two fixed-dimension perturbative schemes: the massive 
zero-momentum scheme and the minimal-subtraction scheme without e ex- 
pansion. The analysis of the corresponding six-loop and five-loop series shows 
the presence of a stable fixed point in the relevant coupling region. Therefore, 
we predict the transition to be continuous. We also compute critical expo- 
nents. The specific-heat exponent a is estimated as a = 0.20(15), while the 
magnetic susceptibility and magnetization exponents 7# and f3u for 3 He are 
1H = -0.34(5), H = 1.07(9). 

PACS Numbers: 67.57.Bc, 64.60.Fr, 05.70.Jk, ll.10.Kk 
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I. INTRODUCTION AND SUMMARY 



The superfluid phases in 3 He have been extensively studied in the years [1-4]. In the 
absence of magnetic field two superfluid phases, named A and B respectively, have been 
identified. However, if one takes into account the weak dipolar interaction [5,1], a new 
phase, called planar phase, should separate the normal state from the A and B superfluid 
states [6,7]. At present, there is no experimental evidence for a normal-to-planar transition. 
This is due to the fact that dipolar interactions are expected to become relevant only within 
the interval \t\ < 2 x 1CT 6 , where t = 1 — T/T c . For pressures between 10 and 30 bar, T c 
varies between 2 and 2.6 mK, so that the interesting temperature interval is approximately 
\T— T c \ PS (4-5)-10 -9 K. Outside this interval normal Fermi liquid or A/B superfluid behavior 
should be observed. Despite the narrowness of such temperature interval, Ref. [6] argued that 
it should be possible to observe the planar phase if the temperature T can be controlled 
to high precision, as it has been done recently for 4 He [8]. It is therefore interesting to 
investigate its nature in detail. 

The normal-to-superfluid transition in 3 He is supposed to be driven by p-wave spin-triplet 
pairing states. Near the transition the critical behavior is controlled by an order parameter 
that is a complex 3x3 matrix, see, e.g., Refs. [9-11,1]. The weak dipolar interaction breaks 
the original SO(3) xSO(3) xU(l) symmetry down to SO(3)xU(l) and the corresponding 
order parameter becomes a complex three-component field ip. If we neglect strain gradient 
effects [12], which are expected to be irrelevant at the transition [6], the Landau-Ginzburg- 
Wilson Hamiltonian describing the critical modes at the normal-to-planar transition can be 
written as [6,7,13] 



/a 



n= / d d x 



(i) 



with /o > and g < 0. Approximate estimates [6,14-16] give go/ fo ~ —1/3. Note that 
stability requires fo > and fo + go > 0. When g < the Hamiltonian is minimized by 
fields ip a = e l6 n a , where n a is a constant real vector, and the corresponding ground-state 
manifold is (S 2 x S ,1 )/Z 2) where S n is the n-dimensional sphere. The case g > with a two- 
component and a three-component field ip is also physically interesting because it describes 
the critical properties of frustrated spin models with noncollinear order [17,18]. The same 
Hamiltonian has also been considered to discuss the critical properties of Mott insulators 
[19]. It should also be relevant for two-dimensional quantum phase transitions in cuprate 
superconductors [20] and for transitions in heavy-fermion compounds such as UPt3 [21]. 

In the mean-field approximation the model with Hamiltonian (1) has a continuous second- 
order transition. On the other hand, e-expansion calculations for an A-component field ip 
find [6,7,17,22] a stable fixed point (FP) in the region g < only for 2-e < N < 2.20-0.57e, 
e = 4 — d. Apparently no FP is found for N = 3, suggesting a fluctuation- induced first-order 
transition. This scenario has been further supported by three-loop calculations within a 
three-dimensional (3-d) massive zero-momentum (MZM) scheme [23] and by a nonpertur- 
bative renormalization-group (RG) study of the effective average action in the lowest-order 
approximation of the derivative expansion [24]. As we shall argue below mentioning a 
few specific physical examples, these results may not be conclusive. For instance, in some 
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physically interesting cases, low-order, and in some cases also high-order, e-expansion cal- 
culations fail to provide the correct physical picture. The location and the stability of the 
FP's may drastically change approaching d — 3, and new FP's, not present for d ~ 4, may 
appear in three dimensions. We mention the Ginzburg-Landau model of superconductors, 
in which a complex scalar field couples to a gauge field. One-loop e-expansion calculations 
[25] indicate that no stable FP exists unless the number N of real components of the scalar 
field is larger than N c = 365 + 0(e). This number is much larger than the physical value 
N = 2. Consequently, a first-order transition was always expected [25]. Later, exploit- 
ing 3-d theoretical approaches (see, e.g., Ref. [26]) and Monte Carlo simulations (see, e.g., 
Ref. [27]), it was realized that 3-d systems described by the Ginzburg-Landau model can 
also undergo a continuous transition — this implies the presence of a stable FP in the 3- 
d Ginzburg-Landau theory — in agreement with experiments [28]. A similar phenomenon 
probably occurs — but this is still a controversial issue, see, e.g., Refs. [18,29,30] — in model 
(1) when g > 0. e-expansion calculations find a stable FP with attraction domain in the 
region g > only for N > N c with [17,22,31,32] N c = 21.80 + 0(e), which seems to exclude 
the physically interesting cases N = 2, 3. These conclusions were apparently confirmed by 
three-loop 3-d perturbative calculations in the MZM scheme [23] and by a nonperturbative 
RG study of the effective average action [33], using approximations based on the lowest 
orders of the derivative expansion. On the other hand, recent studies based on 3-d high- 
order field-theoretical calculations found a stable FP with attraction domain in the region 
go > for both N = 2 and N = 3. Two different renormalization schemes were used: 
the MZM scheme (six loops) [34,35] and the minimal-subtraction (MS) scheme without e 
expansion (five loops) [36]. These results confirm the existence of a new chiral universality 
class and are in agreement with most experiments that observe continuous transitions in 
stacked triangular antiferromagnets [18,29,37]. 

In this paper we investigate the nature of the normal-to-planar transition in 3 He by 
exploiting two different 3-d perturbative field-theoretical approaches based on the Hamil- 
tonian (1). The first one, the MZM scheme, is defined in the massive (disordered) phase 
[38] . It is based on a zero-momentum renormalization procedure and an expansion in powers 
of zero-momentum quartic couplings. As already mentioned, this perturbative scheme was 
considered in Ref. [23], where three-loop series were analyzed without finding evidence for 
stable FP's. The second scheme is the so-called MS scheme without e expansion [39], which 
is defined in the massless critical theory. This scheme is strictly related to the e expansion, 
but, unlike it, no expansion in powers of e is performed. Indeed, e is set to its physical value 
e = 1, and one analyzes the expansions in powers of the renormalized quartic couplings. 

Here we shall present an analysis of the six-loop MZM series that have been computed in 
Ref. [34] and of the five-loop MS series that have been computed in Ref. [32]. We find a stable 
FP in the region g < whose attraction domain includes the region where <?o//o ~ ~~ 1/3, 
which should be the relevant one for the transition in 3 He [16]. The two perturbative schemes 
give consistent results, providing a nontrivial crosscheck of the results. We thus predict a 
continuous normal-to-planar transition, belonging to a new universality class. Our result 
contradicts earlier theoretical works, see, e.g., Refs. [6,23,24]. 

We also determine the standard critical exponents for model (1). We obtain v = 0.59(4), 
v = 0.079(7), 7 = 1.14(8), (3 = 0.32(3) in the MZM scheme and v = 0.63(8), rj = 0.086(24), 
7 = 1.20(15), (3 = 0.34(4) in the MS scheme. A weighted average gives 
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= 0.60(5), 77 = 0.08(1), 7 = 1.16(10), = 0.33(4). (2) 



The errors are such to include the MZM results (with their errors) that apparently are 
the most precise ones. By using hyperscaling we also obtain a = 0.20(15). Note that, 
while v and a are indeed the critical exponents associated with the specific heat and with 
the correlation length, 7, 77, and (3 are not directly accessible in 3 He and are not related 
with the magnetic critical behavior of 3 He. Indeed, the operator that couples with the 
external magnetic field in 3 He is the so-called chiral operator [17,18,40] C a = i^abc^t^c- We 
compute the exponents 7^ and (5h that characterize the singular behavior of the magnetic 
susceptibility and of the magnetization, and r/ H , related to the large- momentum behavior 
of the magnetic structure factor (in field-theoretical terms it is defined in terms of the 
Fourier transform of (C a (0)C a (x)) that behaves as q~ 2+r > H for large momenta q). We find 
r] H = 2.72(7), 7^ = -0.34(5), (3 H = 1.05(7) in the MZM scheme and rj H = 2.63(19), 
7h = —0.34(11), [3h = 1.11(12) in the minimal-subtraction scheme. A weighted average 
gives 

rj H = 2.70(9), lH = -0.34(5), (3 H = 1.07(9), (3) 

where the errors have been computed as before. Note that 7^ is negative so that in the 
high-temperature phase the analytic background gives the dominant contribution to the 
magnetic susceptibility x- Close to T c , we have therefore x ~ const + (T c — T)°- 34 . In the 
low-temperature phase instead x is infinite due to the presence of Goldstone modes. 

The paper is organized as follows. In Sec. II we present the analysis of the six-loop series 
in the MZM scheme. It provides a rather robust evidence for the presence of a stable FP 
with attraction domain in the region g < 0, which is the relevant one for the 3 He superfluid 
transition. This result is fully confirmed by the analysis of the five-loop series in the massless 
MS scheme that is presented in Sec. III. We report the perturbative expansions analyzed 
here in the MZM and MS schemes in App. A and B, respectively. 

II. THE MASSIVE ZERO-MOMENTUM SCHEME 

In this Section we consider the fixed-dimension massive zero-momentum (MZM) scheme 
that describes the disordered massive phase. The relevant RG functions were computed to 
six loops in Refs. [34,40]. In App. A we report the series for the relevant case N = 3. 

The perturbative expansions are asymptotic. Therefore, perturbative series must be 
resummed in order to study the RG flow. Consider a generic quantity S(f, g) = c ijf % 9^ ■ 
In order to determine S(f,g), one may resum the series 

S(xf,xg) = J2s k (f,g)x k (4) 
k 

in powers of x and then evaluate it at x — 1. This can be done by using the conformal- 
mapping method [41,42], which exploits the knowledge of the large-order behavior of the 
expansion, or the so-called Pade-Borel method. The large-order behavior of the coefficients 
is generally given by 
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**(/, 9) ~ g)) k k b [l + OCA:- 1 )] . (5) 

Using semiclassical arguments, one can argue that [34] the expansion is Borel summable for 

/>0, f + g>0. (6) 

In this domain we have 



A(f,g)=aMax[f,f + g]>0, (7) 

where a ~ 0.14777422. Under the additional assumption that all the singularities of the 
Borel transform lie on the negative axis, the conformal-mapping method turns the original 
expansion into a convergent one [41] when conditions (6) are satisfied. Outside this region 
the expansion is not Borel summable. However, if the condition 

f + \g>^ (8) 

holds, then the singularity of the Borel transform that is closest to the origin is still in the 
negative axis. Therefore, the leading large-order behavior is still given by (5) with A(f,g) 
given by (7). In this case, the conformal-mapping method is still able to take into account 
the leading large-order behavior, although it does not provide a convergent series. Therefore, 
one may hope to get an asymptotic expansion with a milder behavior, which still provides 
reliable results. 

The RG flow of the theory is determined by the fixed-points (FP's). Two FP's are easily 
identified: the Gaussian FP (/* = g* = 0) and the 0(6) FP at [43] /* « 0.86, g* = 0. 
The results of Ref. [44] on the stability of the three-dimensional O(M) -symmetric FP's 
under generic perturbations can be used to prove that the 0(6) FP is unstable. Indeed, the 
perturbation term \i/j 2 \ 2 in Hamiltonian (1) is a particular combination of quartic operators 
transforming as the spin-0 and spin-4 representations of the 0(6) group, and any spin-4 
quartic perturbation is relevant at the O(M) FP for M > 3 [44], since its RG dimension y 4;4 
is positive for M > 3. In particular, y AA « 0.27 at the 0(6) FP [44]. 

The analyses of the six-loop series reported in Refs. [34,35] provided a rather robust 
evidence of the presence of a stable FP with g* > and attraction domain in the region 
go > 0, at /* = —0.4(2) and g* = 1.5(1). As already mentioned in the introduction, this 
FP should descrive continuous transitions in frustrated three-component spin models with 
noncollinear order. However, it is of no relevance for the critical behavior of 3 He. Indeed, 
due to the presence of the unstable 0(6) FP, the g — axis plays the role of a separatrix, 
and thus the RG flow corresponding to go < cannot cross the g — axis. Therefore, the 
RG flow relevant for the normal-to-planar transition is determined by FP's in the region 
g < 0. We shall now show that the RG flow obtained by resumming the six-loop /^-functions 
provides a rather robust evidence for another stable FP in the region g < with attraction 
domain including the region go/ fo ~ — 1/3, which should be the relevant one for the 3 He 
superfluid transition. 

In order to investigate the RG flow in the region go < 0, we apply the same analysis 
of Refs. [34,45]. We refer to these references for details. We resum the perturbative series 
by means of the conformal-mapping method [41] that takes into account the large-order 
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FIG. 1. Zeroes of the /3-functions in the f-g plane in the MZM scheme. Pluses and crosses corre- 
spond to zeroes of j3f and f3 g , respectively. The FP coordinates are: /* = 1.08(12), g* = —1.46(10). 
Different colors are used to mark the regions: (1) / + g > 0, (2) / + g < and / + \g > 0, (3) 

/ + \g < o. 

behavior of the perturbative series. To understand the systematic errors, we vary two dif- 
ferent parameters, b and a (see Refs. [34,45] for definitions), in the analysis. We also apply 
this method for those values of / and g for which the series are not Borel summable but 
still satisfy / + g/2 > 0. As already discussed, although in this case the sequence of ap- 
proximations given by the conformal-mapping method is only asymptotic, it should provide 
reasonable estimates, since we are taking into account the leading large-order behavior. In 
Fig. 1 we report our results for the zeroes of the /^-functions, obtained from the analysis of 
the /-loop series, I = 3,4,5,6. For each /3-function we consider 18 different approximants of 
/3(f, g)j f with b = 3, 6, . . . , 18 and a = 0, 2, 4 and we determine the lines in the (/, g) plane 
on which they vanish. Then, we divide the domain < / < 2 and —2.5 < g < into 40 x 50 
rectangles, marking those in which at least three approximants of each (5 function vanish. 
No evidence for a FP appears at three loops, consistently with Ref. [23]. As the number I 
of loops increases, a new stable FP — quite stable with respect to / — clearly appears. Such 
zeroes appear in 9%, 22%, 59%, 89% of the cases we consider for / = 3,4,5,6. Clearly, the 
new zero is increasingly stable as / increases. At 6 loops we obtain a robust evidence for a 



6 



FP at 



f* = 1.08(12), g* = -1.46(10), (9) 

where the error bars correspond to 2 standard deviations. All zeros of the approximants 
with 3 < b < 18 and 2 < a < 4 lie within half of the reported confidence interval and 
only a few with a = give significantly different estimates for /* and g*. Moreover, for 
this smaller class of approximants a stable FP appears in 100% (73%) of the cases for / = 6 
(1 — 5). Estimate (9) is stable with respect to the number of loops. For / = 5 we obtain 
/* = 1.12(18), g* = —1.47(18). Notice that the FP belongs to the region in which the 
series are not Borel summable, but the FP values still satisfy /* + g*/2 > 0. Therefore, 
our resummations should be reliable: the stability of the results with respect to I supports 
this expectation. The FP (9) is stable, since the eigenvalues of the corresponding stability 
matrix have positive real part. Our estimate is ou± = 1.03(29 + 15) ± iO. 72(11 + 10), where 
the first error takes into account the variation of the estimate with a and b, while the second 
one takes into account the error on /* and g*. Note that the eigenvalues are apparently 
complex, suggesting that the FP is actually a focus, as it happens for the chiral FP [35]. 
Using the RG functions rj^ and r) t , we have generated perturbative series with FP value 
corresponding to v, rj, 7, and (5. Correspondingly, we obtain the following estimates of the 
critical exponents: 

v = 0.59(3 + 1), r] = 0.079(2 + 5), 7 = 1.14(6 + 2), (5 = 0.32(2 + 1), (10) 

where, as before, the first error gives the dependence on a and b and the second is related 
to the uncertainty on the FP coordinates. Consistent results are obtained from the analysis 
of the perturbative series for 1/V, r\j f 2 , f 2 /r/, I/7, and 1/(3. The RG relations among the 
different critical exponents are also well verified. 

The exponents reported in (10) follow the usual nomenclature of magnetic systems. How- 
ever, while v is indeed the correlation-length exponent in 3 He, rj, 7, and f3 are not accessible 
in 3 He and in particular, are not related to the magnetic properties of the superfluid. Indeed, 
in 3 He an external magnetic field H gives rise to a coupling of the form [46] 

iJ2 e abcH a 1pl^ c . (11) 

abc 

Therefore, the magnetic exponents are obtained by computing the RG dimension of the 
operator 

C a = ie a bciptipc (12) 

Such an operator was discussed at length in Refs. [17,18,40] in the context of the frustrated 
chiral models and its associated RG dimension was computed at the chiral FP with g > 0. 
In the present case the analysis of the six-loop series of Ref. [40] at the FP (9) gives 

7^ = 2.72(3 + 4), 7^ = -0.34(2 + 3), (3 H = 1.05(4 + 3), (13) 

where the two errors are related to the dependence on a and b and to the uncertainty on 
the FP coordinates. We have also resummed the series corresponding to l/r/ H , 1/ ' f3 H) and 
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FIG. 2. RG flow in the f-g plane in the MZM scheme. It is obtained by using the approximant 
(b = 6, a = 2) for /?/(/, g)/ f and (3 g (f, g)/f- The thick line is related to the limiting RG trajectory 
obtained for go — > The other trajectories reaching the FP correspond to (from top to bottom) 
go /f = -0.085, -0.172, -0.261, -0.354, -0.454, -0.561. Different colors are used to mark the 
regions: (1) / + g > 0, (2) / + g < and / + \g > 0, (3) / + \g < 0. 

1/jh- In the first two cases, the results are in perfect agreement with the estimates (13); 
only for 7^ do we find a significant difference. Note, however, that scaling relations give 
1h = ^(2 — t)h) = —0.42(5) and 7^ = 3u — 2j3 H = —0.33(18), and thus confirm the direct 
estimate (13). 

In Fig. 2 we show the RG flow in the (/, g) plane for several values of the Hamiltonian 
parameters fo and g and for a given approximant of the two (5 functions (other approximants 
give similar results). We consider values in the stability region fo>0 and f + g > and 
restrict the calculation to the region of interest go < 0. We refer to Ref. [47] for details on the 
determination of the RG trajectories and their relation with the quartic parameters /o and 
go of the Hamiltonian. All trajectories corresponding to —0.65 < go/ fo < lie completely 
in the region / + \g > in which resummations should be reliable and are attracted by 
the stable FP. Therefore, the attraction domain of the FP (9) safely includes the region 
9o/fo ~ — 1/3, which should be the relevant one for the normal-to-planar transition in 3 He. 
In Fig. 2 we also report (thick line) the RG trajectory that starts at the 0(6) unstable FP 
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and that can be obtained as the limit go — > 0~ of the trajectories starting at the Gaussian 
FP. Indeed, in this limit the flow first runs close to the g = axis and then follows the thick 
line. This special trajectory is relevant for the analysis of the crossover behavior between 
the 0(6) FP and the stable FP with g < 0; see, e.g., Ref. [47]. 

It is interesting to note that the estimates of the critical exponents (10) are close to 
those associated with the chiral FP at /* = -0.3(1), g* = 1.5(1), i.e. [34] v = 0.55(3) and 
i] = 0.10(1). The reason is that the functions flf, g (f,g) and g) have the approximate 

symmetry 

9^9 = -9, f^f = f + 9, (14) 

i.e. satisfy r)^ t (f,g) fx r)^ >t (f + g,-g), Pf(f,g) ~ Pf(f + g,-g) + P g (f + g,-g), and 
Pg(fi9) ~ —(3g(f + 9, —g)- Perturbatively, the violations of these relations are much smaller 
than the RG functions themselves and appear at three loops in the /^-functions and at five 
loops in T]^ t f Moreover, this symmetry is satisfied by the constant A(f,g), cf. Eq. (7), 
that controls the large-order behavior. This means that the presence of the chiral FP [34] 
/* h fx —0.3, g* h fx 1.5, implies the presence of another FP for /* fx /* h + g* h fx 1.2 and 
g* fx —g* h fx —1.5, in substantial agreement with Eq. (9), with approximately the same 
exponents (10). Note that the approximate symmetry (14) does not hold for r]H(f,g) and 
thus the exponents t]h, Ph, and 7# differ significantly from the corresponding exponents at 
the chiral FP computed in Ref. [40]. 

In conclusion, the analysis of six-loop series in the MZM scheme provides a rather robust 
evidence for a stable FP with attraction domain in the region g < 0, and therefore for the 
existence of a universality class associated with the normal-to-planar superfluid transition in 
3 He. This result contradicts earlier perturbative studies based on e expansion [22] and lower 
order (three-loop) calculations within the same scheme [23], and also a nonperturbative 
study of the RG flow of the effective average action in the lowest order of the derivative 
expansion [24]. 

III. THE MASSLESS MINIMAL-SUBTRACTION SCHEME 

In this perturbative approach one considers the massless critical theory in dimensional 
regularization and renormalizes it in the minimal-subtraction (MS) scheme [48]. In the 
standard e-expansion scheme [49] the FP's, i.e. the common zeroes of the /5-functions, are 
determined perturbatively by expanding in powers of e. Once the expansion of the FP values 
/* and g* is available, exponents are obtained by expanding the RG functions at the FP in 
powers of e. The MS scheme without e expansion [39] is strictly related. The RG functions 
[3f t9 and rj^j are the MS functions. However, e is no longer considered as a small quantity 
but it is set to its physical value, i.e. in three dimensions one simply sets e = 1. Then, 
the FP values /*, g* are determined from the common zeroes of the resummed (3 functions. 
Critical exponents are determined by evaluating the resummed RG functions r/^ and r) t at 
/* and g*. Notice that /* and g* have nothing to do with the FP values obtained within 
the MZM scheme, since / and g indicate different quantities in the two schemes. 

The /9-functions and the RG functions associated with the standard exponents have been 
computed to five loops in Ref. [32] for generic values of N. Moreover, we computed the five- 
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FIG. 3. Zeroes of the /3-functions in the f-g plane in the MS scheme. Pluses and crosses corre- 
spond to zeroes of (3f and (3 g , respectively. The FP coordinates are: /* = 0.89(16), g* = —0.90(24). 
Different colors are used to mark the regions: (1) / + g > 0, (2) / + g < and / + \g > 0, (3) 

/ + h < o. 

loop series of the RG dimension of the operator (12) associated with the magnetic field in 
3 He. In App. B we report the series for the relevant case N = 3. 

The series are resummed by using the conformal-mapping method. Semiclassical argu- 
ments allow us to compute the large-order behavior of the perturbative series. We obtain 
the same formulae reported in Sec. II, cf. Eqs. (5) and (7), with the only difference that 
a = 1/2. We should mention that in this approach the series are essentially four-dimensional. 
Thus, they may be affected by renormalons which make the expansion non-Borel summable 
for any / and g and are not detected by a semiclassical analysis, see, e.g., Ref. [50]. This 
problem should also affect the MS series of the 0(iV)-symmetric theories. However, the good 
agreement between the results of their analysis assuming Borel summability [39] and the es- 
timates obtained by other methods indicates that renormalon effects are either very small or 
absent (note that, as shown in Ref. [51], this may occur in some renormalization schemes). 
For example, the analysis [39] of the five-loop MS series (using the conformal mapping 
method and the semiclassical large-order behavior) gives v = 0.629(5) for the Ising model 
and v = 0.667(5) for the XY model, which are in good agreement with the most precise 
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estimates obtained by lattice techniques, v = 0.63012(16) (Ref. [52]) and v = 0.63020(12) 
(Ref. [53]) for the Ising model, and v = 0.67155(27) (Ref. [54]) for the XY universality 
class. On the basis of these results, we will assume renormalon effects to be negligible in the 
analysis of the two- variable series of the theory at hand. 

As in the MZM scheme, two unstable FP's can be easily identified, the Gaussian FP 
and the 0(6) FP along the /-axis. A stable FP is found in the region g > [36], at 
/* = 0.0(1) and g* = 0.9(1), with critical exponents v = 0.61(6) and r] = 0.09(1), in 
reasonable agreement with the results of the MZM scheme. 

Also in this case there is an approximate symmetry g — > g = —g and /—>/ = / + g, 
which is violated at three loops in the /^-functions and at five loops in the standard critical 
exponents, and which is satisfied by the large-order behavior. Thus, we expect a stable FP 
in the region g < at /* 0.9, g* ~ —0.9, with exponents approximately equal to those of 
the chiral FP at g > 0. 

In order to find the zeroes of the /3-functions, we first resummed the series for the 
functions Pf(f,g)/f and (3 g (f, g) / f ■ For each /5-function we considered several approximants 
corresponding to a and b in the range < a < 4 and 3 < b < 18, as in the MZM case. In 
Fig. 3 we show the zeroes of the /^-functions in the region g < 0. A common zero is clearly 
observed at /* ~ 0.9 and g* ~ —0.9. In order to give an estimate of the FP, we considered 
all independent combinations for the two /5-functions. Most combinations, approximately 
92% and 49% for / = 5 and 4 respectively, have a common zero in the region g < —0.3, 
leading to the estimate 

/* = 0.89(16), g* = -0.90(24), (15) 

where we have again reported two standard deviations as error. As before, essentially all 
zeroes lie in half of the quoted interval — only a few corresponding to approximants with 
a = differ significantly. Estimates (15) are also stable with respect to the number of 
loops and in good agreement with those obtained by using the approximate symmetry. We 
also performed a different analysis, in which we resummed Bf(f,g)/f 2 and B g (f,g)/f 2 , 
cf. Eq. (B4), obtaining consistent results. Notice that in this scheme the stable FP is at 
the boundary of the region in which perturbative expansions are Borel summable, since 
f* + g* ~ 0. This gives us further confidence on the reliability of the final results. 

For the stability eigenvalues, we find that in approximately 33% of the cases the eigenval- 
ues turn out to be real, while in 67% of the cases they are complex. If we average the complex 
ones we obtain iv± = 1.2(4 + 4) ± «0.8(2 + 4); real eigenvalues give instead lo\ = 0.9(8 + 7) 
and U2 = 5(5 + 1). Note that the real results fluctuate wildly, indicating that they are 
probably not reliable. On the other hand, the complex estimate is more stable and is in 
good agreement with the MZM result. 

Critical exponents are obtained as before. For the standard exponents we find 

v = 0.63(3 + 5), 77 = 0.086(13 + 11), 7 = 1.20(5 + 10), (3 = 0.34(2 + 2), (16) 
while for the magnetic (chiral) exponents we find 

7^ = 2.63(5 + 14), lH = -0.34(4 + 7), H = 1.11(4 + 8). (17) 
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Gaussian 0(6) 




0.25 0.5 0.75 1 

f 

FIG. 4. RG flow in the f-g plane in the MS scheme. It is obtained by using the approxi- 
mant (b = 10, a = 1) for j3f{f,g)/f and P g {f,g)/f- The thick line represents the RG trajec- 
tory for <7o — ¥ 1 • The other trajectories reaching the FP correspond to (from top to bottom) 
go /f = -0.085, -0.172, -0.261, -0.354, -0.454, -0.561, -0.679, -0.810. Different colors are 
used to mark the regions: (1) / + g > 0, (2) / + g < and / + \g > 0, (3) / + \g < 0. 

The error is given as a sum of two terms, related respectively to the variation with a 
and b (we used < a < 4 and 3 < b < 18) and to the uncertainty of the FP coordinates. 
Estimates (16) and (17) are in good agreement with the MZM ones: in all cases the difference 
is significantly smaller than the quoted errors, which may be an indication that our error 
estimates are rather conservative. 

In Fig. 4 we show the RG flow in the quartic-coupling /, g plane for several values of fo 
and g in the region fo > 0, fo + go > 0, and g < 0. We report the results for a single 
approximant for both f3 functions. Others give qualitatively similar results. All trajectories 
corresponding to —0.9 ^ go/ fo < lie in the region f + \g > and apparently flow towards 
the stable FP. Therefore, the attraction domain of the FP (15) safely includes the region 
9ol fo ~ — 1/3, which should be relevant for the normal-to-planar transition in 3 He. 
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APPENDIX A: THE SIX-LOOP SERIES OF THE MASSIVE 
ZERO-MOMENTUM SCHEME 

The theory is renormalized by introducing a set of zero-momentum renormalization con- 
ditions involving the two- and four-point correlation functions, i.e. 



ritl(o) 



(Al) 



(SaiSbj + SajSbi) + 9SabSij 



where r^ 2 -* and are respectively the Fourier transform of two- and four-point one-particle 
irreducible correlation functions. In addition, one introduces the function Z t (f,g) that is 
defined by the relation r2' 2) (0) = 5 ai Zf\ where T^ 1 ' 2 ) is the one-particle irreducible two- 
point function with an insertion of ip* ■ ip. The FP's of the theory are given by the common 
zeroes of the Callan-Symanzik /3-functions 



W,g) 



m 



dm 



Pg(f,9) = m 



fo,9o 



dg 

dm 



(A2) 



fo,90 



The RG functions rj^ and r) t associated with the critical exponents are defined by 



d\nm 



fo,90 



Their relation with the critical exponents is given by 

v = %(/* , 9*), v = [2 + g*) - vM*, 91V 1 



(A3) 



(A4) 



The RG functions have been computed to six loops in Ref. [34]. Here we report the series 
for N = 3: 



Pf = ~f + I (14 f + 8 / g + 8 g 2 ) - ^ (109 f 3 + 100 f 2 g + 177 / g 2 + 90 g 3 ) + 

+ (0.792654 f 4 + 1.12206 f g + 2.63235 f g 2 + 2.08529 f g 3 + 0.127631 g 4 ) + 
-(1.02168 f + 1.89478 f 4 g + 5.18772 f g 2 + 5.76831 f 2 g 3 + 2.11992 f g 4 + 
+0.243881 g 5 ) + (1.56686 f + 3.64624 f g + 1 1.4687 f 4 g 2 + 16.3586 f 3 g 3 + 
+ 10.9169 f 2 g 4 + 3.75387 f g 5 + 0.480925 g 6 ) - (2.76292 f + 7.70253 f g + 
+26.9939 f g 2 + 46.7321 f 4 g 3 + 44.5897 f g 4 + 25.4728 f g 5 + 
+7.24362 fg 6 + 0.637996 g 7 ), 



(A5) 
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2 16 

Pg = -9 + ^(2 / + 9) ~ ^7^(!27 f + 127 / g + 6 <? 2 ) + (A6) 

+g (0.926499 f + 1.38975 f 2 g + 0.898551 / # 2 + 0.216645 g 3 ) + 

(1.31884 / 4 + 2.63769 f g + 2.87466 f 2 g 2 + 1.55675 f g 3 + 0.179409 g 4 ) + 
+g (2.10866 f + 5.27164 f 4 g + 8.03232 f g 2 + 6.77566 f 2 g 3 + 2.29411 f g 4 + 
+0.193466 g 5 ) - g (3.9354 f + 11.8062 fg + 22.6695 /V + 25.6626 /V + 
+ 14.855 fg 4 + 3.99177// + 0.406903 g 6 ), 



^ = —(2f 2 + 2fg + 3g 2 )+ (A7) 

+ (0.00379233 / 3 + 0.00568849 f 2 g + 0.00893906 /# 2 + 0.00352145 # 3 ) + 
+ (0.00873136 f 4 + 0.0174627 f g + 0.0395962 / 2 # 2 + 0.0308648 / # 3 + 
+0.00346384 g 4 ) - (0.0027221 f + 0.00680525 f 4 g + 0.0226066 / 3 # 2 + 
+0.0271046 f 2 g 3 + 0.00857277 / g 4 + 0.0000120125 # 5 ) + (0.00954181 / 6 + 
+0.0286254 f g + 0.0902133 / 4 # 2 + 0.132717 / 3 # 3 + 0.0958595 f g 4 + 
+0.034245 fg 5 + 0.00413189 g 6 ), 



Tk = -\(*f + 40) + ^(2/ 2 + 2/^ + 3«? 2 ) + (A8) 

-(0.206554 f 3 + 0.309831 / 2 g + 0.467379 / / + 0.182051 g 3 ) + 
+ (0.194175 / 4 + 0.38835 f°g + 0.930785 f 2 g 2 + 0.73661 / g 3 + 
+0.02992 4 ) - (0.291118 / 5 + 0.727795 / 4 # + 1.88748 fg 2 + 
+2.10343 /V + 0.86341 / g 4 + 0.105499 # 5 ) + (0.444686 / 6 + 
+ 1.33406 / 5 ^ + 4.19137 /V + 6.15931 f 3 g 3 + 4.38435 /V + 
+ 1.527// + 0.176483 g % ). 

The critical exponents associated with a real magnetic external field are determined from 
the RG dimension of the operator C a , cf. Eq. (12). In order to compute it, we define the 
renormalization function Z c (f,g) from the one-particle irreducible two-point function r^ c ' 2 ^ 
with an insertion of the operator C a , i.e. 

r( c - 2 )(0)^ = Z- 1 T ai j, (A9) 

where T ai j is an appropriate tensor normalized so that Z c — 1 at tree level. The RG function 

dlnZ, 



Vc(f,g) 



(A10) 

/o,90 



dlam 

was computed to six- loop in Ref. [40]. For the three-component case one obtains 
^ = -^(/-2/ + ^(3/ 2 -2/^-3/) + 
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-(0.0426744 f - 0.108535 f g - 0.0512203 f g 2 + 0.0777653 g 3 ) 
+ (0.0646127 f - 0.0340279 f g - 0.00236481 f g 2 + 0.0275227 f g 3 + 
-0.0251600 g 4 ) - (0.0713318 f - 0.0868275 f 4 g - 0.0504380 f g 2 + 
+0.00425084 f 2 g 3 - 0.0721756 f g 4 - 0.0120727 g 5 ) + (0.125344 f + 
-0.0394498 f g + 0.0763944 f 4 g 2 + 0.0550844 f g 3 - 0.194565 f 2 g 4 + 
-0.0632578 / g 5 + 0.0152157 g 6 ). 

The critical exponents r/ H , j H , and (3h are obtained by 

rj H = 1 - 2 Vc + 2rj, lH = u{2 - Vh ), f3 H = is(l - rj c + rj), (All) 

where r) c = rj c (f*,g*). 



APPENDIX B: THE FIVE-LOOP SERIES IN THE MINIMAL-SUBTRACTION 

SCHEME 



In the MS scheme one sets 



f = A d ^F(f,g), 
go = A d fi e G(f,g), 



(Bl) 



where the renormalization functions Z^, F, and G are determined from the divergent part 
of the two- and four-point correlation functions computed in dimensional regularization. 
They are normalized so that Z^,(f,g) pa 1, F(f,g) ~ /, and G(f,g) ~ g at tree level. Here 
Ad is a c?-dependent constant given by A d = 2V /2 r(d/2)/3. We also introduce a mass 
renormalization constant Z t (f,g) by requiring ZtT^ 1 ' 2 ^ to be finite when expressed in terms 
of / and g. Here r^ 1,2 ^ is the two-point function with an insertion of ip* ■ tp. Once the 
renormalization constants are determined, one computes the (3 functions from 



W,9) = V 



df 
d/j, 



Pg(f,9) = V 



fo,9o 



dg_ 
d/i 



(B2) 



/o,So 



and the critical RG functions and r) t from 

d log Z^j 



rhi>Af>9) = 



d log /I 



(B3) 



Their relation with the critical exponents is given in Eq. (A4). The /^-functions have a 
simple dependence on d: 



P f (f, g) = {d- 4)/ + B f (f, g), /3 g {f, g) = (d- A)g + B g (f, g), 



(B4) 



where the functions Bf(f, g) and B g (f, g) — as well as the RG functions g) — are inde- 

pendent of d. 
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In the following we report the five-loop series [32] for N = 3: 

B f = \ (7 f 2 + 4 / g + 4 g 2 ) - ± (24 f + 22 / 2 g + 39 / / + 20 /) + (B5) 

+ (9.07446 f 4 + 12.5703 / 3 g + 29.4275 / 2 / + 23.8424 f g 3 + 2.41477 /) + 
-(46.7683 f + 85.6528 f 4 g + 227.316 / 3 / + 256.217 / 2 / + 110.777 f g 4 + 
+16.3103 g 5 ) + (296.166 f + 677.379 f g + 2026.83 / 4 / + 2898.93 / 3 / + 
+2096.65 f g 4 + 791.170 / / + 109.827 /), 

5 9 = £ [ (2 / + <?)- ^(28 / 2 + 28 fg + g 2 ) + (11.1573 / 3 + 16.7359 f 2 g + (B6) 

+9.75736 / / + 2.06019 g 3 ) - (62.5357 / 4 + 125.071 f g + 123.278 / 2 / + 
+60.7252 / g 3 + 10.5496 /) + (422.235 / 5 + 1055.59 f 4 g + 1458.99 / 3 / + 
+ 1132.71 / 2 / + 411.714 fg 4 + 52.1814 /) ], 



% = 1 (2 / 2 + 2 / g + 3 /) - ± (14 / 3 + 21 f 2 g + 33 / g 2 + 13 /) + (B7) 

+ 25^2 (86 /4 + 172 /S * + 474 /2 ^ + 388 f 9 " ~ 9 ^ + 
-(0.709208 f + 1.77302 f 4 g + 4.78026 / 3 / + 5.39737 / 2 / + 

2.31434// + 0.325646/), 
Vt = -|(2 / + <?) + ^(2 f 2 + 2 / + 3 g 2 ) + (B8) 

-- 1 -(208 / 3 + 312 / 2 g + 486 / / + 191 g 3 ) + (7.00217 f 4 + 14.0043 / 3 <? + 
lUo 

+29.2911 f g 2 + 22.2890 / g 3 + 3.64149 /) - (33.6900 / 5 + 84.2251 f 4 g + 
+207.371 / 3 / + 226.832 f 2 g 3 + 108.339 / g 4 + 19.1618 /). 

We also computed the MS expansion of the RG dimension of the operator C a , cf. Eq. (12). 
For this purpose, we computed the renormalization constant Z c (f, g) by requiring Z C T^ C ' 2 ^ 
to be finite when expressed in terms of / and g, where r( c ' 2 ) is the one-particle irreducible 
two-point function with an insertion of the operator C a . Then, one defines the corresponding 
RG function 



d\ogZ c 

Vc{f,g) = 



(B9) 

fo,go 



d log// 

The resulting series for the three-component case is given by 

Vc = -\{f -2g) + l(3f 2 -2fg-3g 2 )- ^(22 f 3 -42f 2 g- 27 / g 2 + 26/) + 
3 9 216 

+ (1.93436 f 4 - 2.50866 fg- 2.13390 f 2 g 2 + 0.910908 f g 3 + 0.285515 g 4 ) + 
-(8.73201 f - 9.61852 f 4 g - 9.00230 f 3 g 2 + 1.92888 f 2 g 3 - 0.333945 / g 4 + 
-0.810730/). (B10) 
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The critical exponents rjn, 1h and /3h are computed by using relations (All). The expansion 
of 7] c was computed for generic values of N. We verified that r) c (f, 0) coincides with the 
expansion of the RG function t]t, related to the insertion of the spin-2 operator in the two- 
point function, in the 0(2iV) theory obtained by setting g = (perturbative five- loop series 
for 4> T = (2 + r] T — r])is were computed in Ref. [44]). Moreover, for iV = 2 one can prove and 
check the relation i] c (2x/3, —x/3) = r) t (2x/3, —x/3), and r) ti <t>(2x/3, —x/3) = r]*J where rj*J 
are the RG functions of the XY (0(2)-symmetric) theory, which can be found in Ref. [55]. 
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